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On curve and surface 
stretching in turbulent flow 

By N. ETEMADI 1 

Cocke (1969) proved that in incompressible, isotropic turbulence the average 
material line (material surface) elements increase in comparison with their initial 
values. We obtain rigorously, among other things, good estimates of how much 
they increase in terms of the eigenvalues of the Green deformation tensor. 

Introduction 

In the following note we will study deformation of material curves and surfaces 
convected in a turbulent flow. We will do this by looking at the notion of material 
line and surface elements as a means of generating hypotheses for the so called 
flow path ( Lagrangian ) of the motion. Then we will use these hypotheses to 
obtain upper and lower bounds for the evolution of the ensemble average of the 
arc length (surface area) of an arbitrary curve (surface) in time. 

For the definition of material line and surface elements and their historical 
background see Monin and Yaglom (1975). Cocke (1969) is the first who gen- 
erated these mathematical assumptions, see section 1.1, “implicitly” and gave 
a convincing proof of them for isotropic turbulence. Orszag (1970, 1977) takes 
these assumptions for granted, and by a variant of Cocke’s arguments obtains 
somewhat weaker results than Cocke’s, see Remark 1.2. Our work complements 
the work of Cocke. Namely, we will bring out these assumptions in section 1, 
and we will show, first, that once one accepts these assumptions then Cocke’s 
(1969) results can be improved to obtain “tight” upper and lower bounds for the 
ensemble average of material lines and surfaces. Then, in the remaining part of 
section 1 and in section 2 we carry on the results to arbitrary curves and surfaces, 
and in turn we also obtain upper bounds for moments of the dispersion between 
two points moving in the flow at any given time in terms of their separation at 
initial time. 

Throughout our work we will use x(a, t) as the Lagrangian representation of 
the flow, i.e. the trajectory followed by the particle which is at a at initial time 
to. We will assume that x is smooth enough in a space-time region and is, for 
fixed t, an invertible mapping so that our manipulations are legitimate. We 
will also use |t>| as the magnitude of an arbitrary vector t;, whose components, 
without danger of confusion, will be denoted by Vi,t/2, and t> 3 . 
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1. Curve Stretching 

The motion of an infinitesimal material line 61 is governed by the equation 


d6l d6l ^ c , 

— = («-V)u or — = 


(i) 


with 61 (t 0 ) = S°l. Where ^ is the material derivative following the motion, see 
e.g. Monin and Yaglom (1975) section 24.5. It is easy to check that the above 
equation has the solution, 

3 

61 = ^2x <k 6°l k ;x(o,t 0 ) = a , (xij(a,t 0 )) = (fty) = I, (2) 

k= 1 


where I is the identity matrix. 

Let p(a, t) be the Lagrangian density. It is easy to show that 


det(nj(a, t )) 


p(o» to) £o 
p(o,t) ~ p t ’ 


( 3 ) 


where the left hand side is the determinant of the Jacobian matrix (or the 
Jacobian) of the transformation y = x(a, t),t > to, see Batchelor (1977) p.79. 
In particular, the Jacobian is one when the flow is incompressible. 

Next consider 


\6l\ 2 = 6 0 l T (x i<j ) T (x iJ )6°l = 6 0 l T W6°l. (4) 

Clearly W is a symmetric non-negative definite matrix. Since its determinant is 
(po/pt) 2 , consequently its eigenvalues, say ivi,W 2 ,vJs, are strictly positive and 
we have 

2 

det(W) = det(Xi,j) 2 = iv x w 2 w 3 = (5) 

Pt 

Let A = [a,]) be the unitary (rotation) matrix corresponding to diagonalization 
of W. Thus |i45°/| = |5°/|. We can rewrite, 

\6l\ 2 ( A6°l)f (. A6°l)l ( A6°l)% 

\6°l\ 2 “ Wl \A6°l \ 2 + W2 \A6°l \ 2 + W * \A6°l\ 2 ( 6 ) 

= (sin 2 6 cos 2 ip) w i + (sin 2 Osin 2 ip) + (cos 2 0)ws, 

4 fOl 

where 0, ip are the usual spherical coordinates of the unit vector 

Now we are in a position to take ensemble averages of both sides of this equa- 
tion. To do so, we need to make some assumptions about the joint probability 
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distribution of the variables involved; namely Wi’s, 0 and t p. This is where the 
physics of the problem come in. The following three assumptions have been 
extracted from the work of Cocke (1969) who proved them to be true for incom- 
pressible, isotropic turbulence. 

Assumptions; 

(1) Wi’s are independent of 0 and ip, 

(2) 0 and ip are uniformly distributed over the unit sphere, 

(3) Wi’s are identically distributed. 

Assumption one implies that the random orientation of the element 61 is un- 
correlated with its random deformations along the principal axis. Assumption 
two simply says that 61 is equally likely to be oriented in any direction. Finally, 
assumption three means that the random deformations of 61 along its principal 
axis have the same probability law. 

Proposition 1.1. Under Assumptions (l) and (2) we have, 

^(v^ + V^+V^) < (|^||) < ^<vW+ v/^+V^). (0 

!<(«>! + + «*)*} < (j^jj) < <(u>l + + U/s)*>. («) 


Proof. Since y/x is concave and the coefficients of Wi's add up to one, from 
(6) and (19) we obtain, 

l«l 


i«°/i 

Also it is clear from (6) that, 

i«i 


> (sin^dcos^xJj)y/uy[ + {sin^0sin^rp)y/w2 + [cos 2 6) y /ws. 


i«°'i 


< \sinOcosip\y/wl -I- \sin0sinxp\y/w2 + \cos0\y/w3. 


( 7 ) 


( 8 ) 


Next, use the joint density function of 0 and ip, i.e. j^sinO 0 < 0 < it, 0 < 
ip < 27T, to conclude that the average of the coefficients of tuj’s and their square 
roots, in (6), are | and respectively, and this will in turn give us (*). Since 
the coefficient of to, ’s in (6) are bounded by one we are only left to show the left 
hand side of (tt). To see this note that again by the fact that yfx is concave, 
from (6) and (19) we have, 

| 5 /| , . j.,(sin 2 0cos 2 tp)w 1 + (sin 2 Osin 2 ip) w-2 + (cos 2 0)w3, J 

iToTT = ( w i + w 2 + u>s) 2 p - — ; 

^ , . i . Ism0cos^i|iyi + \sin 0 sinip\w 2 + Icos^Ium. 

> [wi +W2 + t«3)>[ J 1 1 '■ ! ]. 


1* 


lt>l + lt>2 + U>3 


( 9 ) 
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Now the average of the coefficient of Wi's in the above line is | and we are 
through by Assumption (1)./// 

Proposition 1.2. Under Assumptions (1) and (2) and/or (l) and (3) we have, 


j( w i + u>f + w$) < ((i^jj) P } < + wj) fo 0 < p < 2, (t) 


= ^(u>i +™2 + W 3 ) for p = 2, 


(*) 


< ( (M r) < { ^1±A±A } to , < „ < OC, (m) 


^(log{wixV2W 3 )) < 


(iv) 


Proof. In light of what we have said in the proof of the above proposition, 
(u) is immediate under Assumptions (1) and (2). Now we can conclude («) by 
utilizing Assumptions (1) and (3) in (6) to obtain, 

( |j^j| 2 ) = (sin 2 0 cos 2 ip) (tv i) + (sin 2 9 sin 2 tp) (tv 2 ) + (cos 2 0)(w 3 ). 

= (sin 2 0co8 2 tp + sin 2 dsin 2 tp + cos 2 0)(tv i) = (u>i) (10) 

= + tt>2 + W 3 ). 

The left hand side of (t) and the right hand side of (m) follow from an inequal- 
ity like the one given in (7). y/x should be replaced by , and the direction 
of the inequality should be reversed due to convexity when p > 2. The right 
hand side of (t) is trivially true and the left hand side of (m) is a consequence 
of < ((|^r|j) p )» P > 1, see (18), and the right hand side of (*) for p = 1. 

Finally to obtain (iv) use the concave function log(x) rather than yjx in (7) 
under Assumption (l) and (2), and an argument similar to the one given at (10) 
under Assumption (1) and (3 )./// 

Remark 1.1. Note that the left hand side inequalities in the above propositions 
are strict unless W! = W 2 = W 3 . For all the concave functions involved are strictly 
increasing. Clearly this happens only when W is the identity matrix, meaning 
pure rotation. Furthermore, this has to be the case with probability one in order 
to have equality in the above propositions, which is certainly not an interesting 
case. 
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) > 1. 


for any p > 0. 

Proof. Cocke (1969) gives an a priori proof to the effect that the Assumption 
(1) to (3) are true for isotropic incompressible flow. This invokes the above 
proposition and together with Remark 1.1 and the fact that detfaj) = 1, see 
(3), we obtain the first inequality. The second one follows from (18), and the 
one we have just established as follows; 

<1^) = - exp W° 9 ]iwjF>i = ex p{p( l °9 > i- (ii) 

We could have also use d the arithmetic geometric mean inequality, (w[ + + 

u>5)/3 > {/(u>ittf2U>3) r , to achieve the same end./// 

Remark 1.2. The proof of the Assumptions (l) to (3) is the thrust of the work 
in Cocke (1969) in which he has also shown Proposition 1.2 (ii) and (iv), the 
above remark and corollary by using the same argument. Orszag (1970) takes 
these assumptions for granted, follows Cocke’s argument and obtains a weaker 
result, ( ) > 1, see the above corollary for p = 2. Note that this result does 
not imply that the average material line stretches. Now if we agree on all three 
assumptions, then it is trivial to show that (W) = 7(^»y); where 7 is the right 
hand side of (ii) in Proposition 1.2. In this connection see also Orszag (1977), 
p.240-241. 

Next we extend the above results to an arc length following the flow. The 
statement of the inequalities needed to carry this on can be found in the ap- 
pendix. Let C(s;to) : [a, 6] — ► R 3 be a parametric representation of a non- 
random curve at time to. Then C(s;f) = x(C(s;t 0 ),t) is the corresponding 
random curve , following the flow at time t. Let C to and C t , t > to be their arc 
lengths, respectively. For homogenous turbulence define, 


o(() = ( vW+^£+v^ ) 5 m - T 

< 'lit) = \f( 


W1+W2 + U >3 


) 


u>i + u;2 + W3 


) 


where the first and second inequalities are the consequence of (19) and (18) 
respectively. 

Theorem 1.1. For an isotropic, incompressible flow, 

Ct 0 < exp{(log( C t ))} < (C t ) < i(t)C to , 


(0 
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C t0 < a(t)C to < (C t > < ^a(«)C t0 , (it) 

^J-0(t)C to < (C t > < V3/?(t)C t0 . (Hi) 

Proof. The above corollary implies that for the non-random vector 6 °l = 

C"(s;« 0 ), 


(Ml 2>,*C£(s;t 0 )|)> = (Ml«l)> > log(\ 6 °l\) = M|C'(s;t 0 )|). (12) 

fc = l 


Consequently, 


exp{(/og(C t ))> = exp{(log[ j <*«]» 

(by (20), Si = [a, 6]) > J exp{(log (\ 9 - —^’ -|))} ds 

rb $ 

= exp{(log(\Y^ x ,k(C(s;t 0 ),t)C' k (s;t 0 )\)}} ds (13) 

Ja fc= l 

(by (12)) > [ exp{/op(|C'(s;t 0 )|)}ds 

J a 

= [ b \C'(s]t 0 )\ ds = C t0 . 

J a 


Next, inequality in (t) is an immediate consequence of (18) with <j> = exp(x). 
The following one is true by virtue of (18) with <f> = y/x, and Proposition 1.2 
(ii). For the rest of the inequalities all we need is Proposition 1.1./// 

Remark 1.3. Since in an inviscid flow vortex lines remain vortex lines, the 
above theorem is an statement about their evolution in time when the flow is 
isotropic and incompressible. The same is also true for vortex sheets which will 
follow from the discussion in section 2. 

Remark 1.4. The reason for presenting various upper and lower bounds is for 
their potential in applications. For instance, with regard to the above remark and 
under the same conditions, one can compute 7 2 (t) as the ratio of the enstrophy 
at time t to time t 0 . 

Proposition 1.2 will give us “tight” upper and lower bounds for | | p ds) 

in an obvious way. This can be used partially to get information about the mo- 
ments of C t . 
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Theorem 1.2. For an isotropic, incompressible flow, let a p {t) = (wf + w£ + 


ty|)/3. Then, for p > 0, 

(C t0 ) p < <(C t ) p > 

(0 

for 0 < p < 1, 


a p (t)[mm (£[0 ,i]{|C , (s;to)r 1 >] C to < <(C t ) p > < (C t ) p < c p 

(Cto) p (tt) 

for p > 1, 


(C,.)» < c? (C..)' < (C,)" < ((C.)') < 3o,(()(max, < |0, 1 |{|C'( 51 t 0 )|'’-‘}] C«„ 

(*0 


where C 2 (ci) is the minimum (maximum) of the coefficients of C to in the upper 
(lower) bounds for (C t ) in Theorem 1.1, and without loss of generality, we have 
assumed [o,6] = [0,1]. 


Proof. The easiest way to handle (») is to raise the left hand side of (t) in 
Theorem 1.1 to power p, take p inside the log , and then use (18) with <j> = exp(x). 
The left hand side of («) and the right hand side of (m) are also the cosequences 
of (18) with 4> = x p in an obvious way, and the remaining parts are easily followed 
by Theorem 1.1./// 

The next corollary will give us information about the moments of evolution 
of a straight line segment in turbulent flow and also the moments of dispersion 
of two points as time goes on. 

Corollary 1.2. Let di = mm[c5,Op(l)] , ^2 = max[c^,Za p {t)\ with ci,C 2 and 
a p (t) as in the above theorem.Let C(s;$o) = oi + s(a 2 — ai), se[0, l] . Then for 
an isotropic, incompressible turbulence, 

dx|o 2 - ai| p < ((C t ) p ) < d 2 |a 2 - Oi| p for p > 0, (t) 

(|z(o 2 ,f) - x(oi,f)| p ) < d 2 |a 2 - ai| p for p > 0. («) 

Proof. Note that | | = I 02 ~ <*i| and use the above theorem./// 

2. Surface Stretching 

In this section we will extend the above results to the evolution of a surface 
area in turbulent flow. It turns out, just as in the work of Cocke (1969), that 
only minor modifications are needed to do so. 

Let 6° l and 6°k be two infinitesimal material line at t — t 0 . We can form an 
infinitesimal material surface by taking the vector product of these vectors, i.e. 
6°S = 6°l X 6°k. This at time t becomes, 

3 3 

6S = Six 6k = (%2XjS 0 lj) x Xj6°kj). 

3 = 1 3 = 1 


( 14 ) 
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Let the matrices W, A and the eigenvalues, Wi’s, of W, be as before. From the 
identity, 

|t>i x v 2 \ 2 = |vi| 2 |« 2 | 2 - («i * w 2 ) 2 , (15) 

for any given vectors Vi and t >2 and the fact that |6°5| 2 = \A6°l x A6°k\ 2 , A 
being unitary, we can easily obtain, 

\6S\ 2 (A6°l x A6°k)l (A6°l x A6°k )% ( A6°l x A6°k)% 

|fi°S | 2 " W2Ws \AS°l x A6°fc| 2 + WlW 3 I A6°l x A6°k \ 2 + WlW 2 \A6°l x A6°k | 2 

= (sm 2 0cos 2 0)u>2t«3 + (sin 2 Osin 2 tp)wiWz + (cos 2 6)wyW 2 . 

( 16 ) 

Now (16) plays the role of (6), and we only need to modify Assumption (3) as 
follows; 

A ssumption; 

(4) wiw 2 , W 1 W 3 and w 2 w$ are identically distributed. ( Note that Assumptions 
(3) and ( 4 ) are equivalent for incompressible flows.) 

Now all we need to do is to replace wi,w 2 , by w 2 ws,w 5 wi,wiw 2 ; Sl,6°l 
by 6S,6°S , and C t ,C to by S t ,S to , and [a, 6] by D, respectively, in the above 
results, including the remarks, to obtain the new ones corresponding to surfaces. 
Where we let S(u, v; to) be a parametric representation of a nonrandom surface 
on a region D on the plain, S to be its area, and, 


'<-!L 


,aj(g(u,t>;t 0 ),t) dz(S(u,t;;to),t) , 


du 


dv 


dudv , 


(17) 


the area at time t. 

The only nonsymbolic modification of the proofs are: (a) In Theorem 1.2 the 
area of D has to be one or otherwise the right (left) hand side of (m) ((**)) has 
to be multiplied by that area to the power p — 1 due to the correct usage of 
Jensen’s inequality, (b) In Corollary 1.2 the notion of a distance between two 
points has to be replaced by an area of a region on a plane and its left hand side 
of (ii) to be interpreted correctly. 

Remark 2.1 We could have always used Proposition 1.2 and its counterpart 
for surfaces to obtain upper and lower bounds for moments of material lines, 
material surfaces, arc lengths, and surface areas at the expense of having different 
constants. Compare Proposition 1.1 with Proposition 1.2 when p = 1. The 
difference between these two propositions becomes significant if one can realize 
physically non-isotropic incompressible flows that satisfy only one pair of the 
assumptions, involved in these propositions, rather than all of them. Finally, 
for analogous results concerning homogeneous turbulence we invite the reader 
to consult Corrsin (1972). 
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Appendix 

(Jensen's Inequality) Let X be a random variable and (j> a convex (concave) 
function containing the range of X. Assume both X and 4>{X) have ensemble 
averages, then 

•KW) < (>) <*(*)>• (is) 

Proof. See any standard graduate textbook in probability theory or measure 
theory e.g. Billingsley (1986) p.283. 

The following special case of Jensen’s inequality has been used frequently; let 
Pi, P2, and P 3 be three positive numbers whose total sum is one, let aj., 02 , and 
03 be any real numbers. Then with <j> as above we have, 


3 3 

*£>«,) < (>) J^Pi^(a,). 

1=1 1=1 


(19) 


Proof. Let X in (18) be the random variable which takes the value a,- with 
probability p,-, t = 1, 2, 3./ / / 

(Dunford and Schwartz[ 5], p.535). Let (5,£,p) and (5i,Ej,pi) be positive 
measure spaces. Assume p(5) = 1 . Then if K is a p. x pi - measurable function 
defined on S x Si, 


/ exp{ I log\K(s,s 1 )\p(ds)}p 1 (dsi) < exp{ / \log{ / |/ir(s,si)|pi(dsi))]p(<fs)}. 

JSi Js Js J Si 

( 20 ) 
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